In-plane magnetic penetration depth in NbS 2 
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We report on the temperature dependence of the in plane magnetic penetration depth (X a b) 
and first penetration field (Hf oc l/X 2 ab (T) for H\\c) in 2H-NbS2 single crystals. An exponential 
temperature dependence is clearly observed in X a b(T) at low temperature, signing the presence of 
a fully open superconducting gap. This compound is the only superconducting 2H-dichalcogenide 
which does not develop a charge density wave (CDW). However as previously observed in 2H-NbSe2, 
this gap ( Ai = LlksTc) is significantly smaller than the standard BCS weak coupling value. At 
higher temperature, a larger gap (A2 = 1.8 fcsTc) has to be introduced to describe the data which 
are compatible with a two gap model. The superconducting gaps are hence very similar in NbS2 
and NbSe2 and we show here that both of them open in the strongly coupled Nb tubular sheets 
independently of the presence of a CDW or not. 



I. INTRODUCTION 

In numerous strongly correlated systems, superconduc- 
tivity appears in the vicinity of a second electronic in- 
stability. These superconductors then often show devi- 
ations from the weak coupling BCS theory which raises 
the question of the interplay between the two orders^. 
This question is particularly well illustrated in dichalco- 
genides. Indeed, in lT-TiSe^, 2H-TaS 2 45 , or 2H- 
NbSe2 superconductivity coexists with a charge density 
wave (CDW). In 2H-NbSe2 the charge density wave de- 
velops for instance at Tjj = 33 K and coexists with su- 
perconductivity (SC) below T c = 7.1 K. 

In parallele, specific heat, tunneling spectroscopy as 
well as de Haas van Alphen measurements^"— showed 
that NbSe2 is characterized by the coexistence of two su- 
perconducting gaps : Ai = 1.0 ± 0.1/cbT c which is much 
smaller than the BCS weak coupling value (A//cbT c ~ 
1.7) and A2 = 2.0 ± 0.3fc#2c- It was then initially sug- 
gested that the CDW might be at the origin of Ai— . 

However, scanning tunneling spectroscopy (STS) and 
specific heat measurements clearly indicated the pres- 
ence of a reduced superconducting gap (of similar am- 
plitude) in the NbS2 system as well 15,16 showing that 
this interpretation does not hold the ground. Indeed, 
this iso-electronic compound is particularly interesting 
to probe the possible influence of the CDW on the 
superconducting states in dichalcogenides, as it is the 
only system among the four 2H layered superconduct- 
ing transition-metal dichalcogenides 17 MX2 (M=Nb,Ta; 
X=S, Se) which does not exhibit a CDW transition. 
Those measurements hence clearly showed that the pres- 
ence of two superconducting gaps in 2H-dichalcogenides 
is not a direct consequence of the presence of a CDW in 
some of them. 

In a first scenario, those two gaps can be considered 
as being the limiting values of one single anisotropic su- 
perconducting gap (without nodes) whereas in an alter- 
native scenario they can be related to different values of 
the electron phonon coupling constants (X e - P h) on differ- 



ent sheets of the Fermi surface (FS) (so called multiband 
superconductivit}^) . 

Clear evidence for multiband superconductivity were 
first obtained in MgB2 for which X e - P h varies from ~ 1.2 
for the 2D-a band to - 0.4 for the 3D-tt bandit. Since 
the FS of NbSe2 is composed of four strongly coupled 
quasi-2D Nb-derived cylinders (bands 17 and 18 with 
X e - P h ranging from 0.8 to 1.9) and one Se-derived pan- 
cake sheet with a much smaller X e - P h value ~ 0.3 (band 
16) 8 , it was tempting to ascribe the small superconduct- 
ing gap to this reduced X e - P h 21}22 However ARPES 
measurement s 1 1 1 14 suggested that, in NbSe2, both gaps 
open on the tubular sheets and the situation iseems to 
be quite different from that of MgB2 as the two gaps 
could both be associated to the same (strongly coupled) 
FS sheets. However, an experimental confirmation that 
the tow gaps are also related to the same FS sheets in 
NbS2 was still lacking. 

In order to address this issue, we have performed in 
plane penetration depth and penetration field measure- 
ments in NbS2 single crystals. Indeed, A is a directional 
probe and corresponding measurements are only sensitive 
to supercurrents flowing perpendicularly to the direction 
of the applied magnetic field. Band structure calculations 
in NbSe2 have for instance shown that the tubular sheets 
contribute to 98% of the transport current in the basal 
planed, and the present measurements will hence only 
probe the superconducting properties of those strongly 
couple sheets. We show here that the temperature de- 
pendence of X a b are remarkably identical in NbSe2 and 
NbS2- This temperature dependence clearly confirms the 
presence of a small superconducting gap Ai = 1.1 fee T c 
(being much smaller than the standard BCS weak cou- 
pling value) and that, as previously observed in NbSe2, 
this gap opens on as least one of the Nb sheets. At higher 
temperature the superfluid density can be deduced from 
A(T) introducing a calibration factor (see below) which 
can be accurately obtained combining A(T) and first pen- 
etration field measurements. We then also show that a 
second gap value is required to describe the tempera- 



ture dependence of this superfluid density which can be 
well described by a two gap model introducing gap values 
very similar to those previously obtained in specific heat 
measurements. We hence confirm that the presence of 
two well defined superconducting gap values, both open- 
ing on the Nb tubular sheets in a generic feature of 2H- 
dichalconenides. 

STS and specific heat (C p ) measurements clearly 
showed that NbSe2 and NbS2 systems share the par- 
ticularity of having both two superconducting gaps but 
also highlighted some fundamental differences. Indeed, 
whereas vortices exhibit a clear stellar shape in NbSe2, 
they appear to be isotropic in NbS2- Moreover, the two 
systems have very different anisotropics of their upper 
critical fields increasing from ~ 3.2 in NbSe2 9 to ~ 8.1 in 
NbS2 24 . However, the vortex shape and H C 2 values are 
determined by the coherence length (£) which is related 
to the superconducting gap but also to electronic prop- 
erties of the normal state (Fermi velocity and possibly 
electronic mean free path in the so-called dirty limit). 
On the contrary the temperature dependence of the pen- 
etration depth is determined only by the superconducting 
gap structure. Those later measurements hence provide 
very valuable information to understand those differences 
which will be discussed at the end of the paper. 



II. SAMPLE AND EXPERIMENTS 

2H-NbS2 samples were grown using the vapour trans- 
port growth technique with a large sulphur excess, as 
described elsewhere 25 . Single crystals have a flat hexag- 
onal or triangular shape and typical dimensions of 150 2 x 
30 jam. The large flat faces grow perpendicular to the c- 
axis. The 2H crystallographic structure has been checked 
on each sample using a four-circle diffractometer. The 
bulk critical temperature was checked on several sam- 
ples of each batch by ac specific heat 16 . All the mea- 
sured crystals showed sharp superconducting transitions 
(AT C ~ 0.3 K) at T c = 6.05±0.2 K emphasizing the excel- 
lent bulk homogeneity of each crystal. Previous scanning 
tunnelling microscopy measurements on single crystals of 
the same origin confirmed the absence of CDW down to 
0.1 K— . To avoid a contamination by oxygen intercala- 
tion between layers which would reduce T c and smears 
the transition 24 , samples were kept in vacuum. All the 
measurements were carried out in the first month follow- 
ing the growth of the samples. 

A high stability LC oscillator operating at 14 MHz, 
driven by a Tunnel Diode has been used to measure 
the variation of the magnetic penetration depth. The 
tiny sample, placed in the center of the solenoid coil, acts 
as a small perturbation to the inductance (L). By chang- 
ing the temperature of the sample, the magnetic penetra- 
tion depth variations produce a change of the inductance, 
leading to a change in the LC resonant frequency Sf(T). 
The amplitude of the AC excitation field is smaller than 
1/aT (that is «H c i) to keep the sample in the Meissner 
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Figure 1: Temperature dependence of the frequency shift 
5f (T), for B (square) D (triangle) (right scale) and C8 (cir- 
cle) (left scale) NbS2 single crystals. Inset: microscope images 
of each sample (note the bigger scale for C8). 

state. This technique enables to measure the variation 
of A with an excellent sensitivity as our setup can detect 
relative frequency changes, Sf/fo of a few 10 -9 , which 
corresponds to changes in A on the order of 1 Ain mm 3 
samples. The sample is glued with vacuum grease at 
the bottom of a sapphire cold finger whose temperature 
is regulated between T m ^ n =0.5 K and 10 K whereas the 
LC oscillator remains at fixed temperature (3K). The fre- 
quency shift Sf=f(T) — f(Tmi n ) is then proportional to 
the magnetic penetration depth : 

S\(T) = R 5 -^j± (1) 

where A/o is the change induced by the introduction of 
the superconducting sample in the coil (measured in-situ 
by extracting the cold finger from the coil at the low- 
est temperature), and R an effective dimension which 
depends on the sample geometry. In principle, R can 
be evaluated by using the formula proposed by Prozorov 
and al. 29 in the case of cylinders. However, recent mea- 
surements in pnictides highlighted the difficulties of eval- 
uating this effective dimension, in particular for mag- 
netic fields perpendicular to the ab-planes (main sample 
surface)^. Indeed, very different SX(T) were obtained by 
different groups, most probably due to the roughness of 
the sample edges which undermines this evaluation. 

To overcome this difficulty we have performed first pen- 
etration field (Hf) measurements using miniature GaAs- 
based quantum well Hall probe arrays. The remanent 
local DC field (B rem (H a )) in the sample has been mea- 
sured after applying a magnetic field H a and sweeping 
the field back to zero. In the Meissner state, no vortices 
penetrate the sample and B rem remains equal to zero up 
to H a = Hf. A finite B rem value is then obtained for 
H a > Hf due to vortices which remain pinned in the 
sample. As Hf is directly proportional to the first criti- 
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cal field H c \ oc 1/ A 2 , the normalized in plane superfluid 
density : 

A 2 1 
p ab (T) = ^JL_ = _ _ _n_ §f{T)/)2 (2) 

is expected to vary as Hf(T)/Hf(0) (Ao being the zero 
temperature penetration depth). Even though the un- 
certainties on the Hf values hinders any precise determi- 
nation of the gap values, those measurements can then 
be used as milestones to obtain the R/XoAfo coefficient 
(see inset of Fig. 3). 

III. EVIDENCE FOR A SMALL GAP VALUE IN 

NBS 2 
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Fig. [H shows typical Sf(T) curves between 0.08 T c 
and 0.85 T c for three NbS2 single crystals (see inset of 
Fig. [1]). The magnetic field is applied perpendicularly to 
the layers, i.e. along the crystallographic c axis, so that 
super currents are flowing in the layers, probing the in- 
plane penetration depth, X a ^. Samples B and D were 
grown in the same batch, and have, in-plane dimensions 
equal to 110 pm and 140/im, respectively and a thick- 
ness t = 30/im (resp. 25/xm). Sample C8 was grown in 
a second batch and has a larger size (430 2 x 110 jam 3 ). 
The frequency shift is ten times higher for C8 than for 
other samples due to it larger volume but all renormal- 
ized curves are fully identical underlying the excellent 
reproducibility of the temperature dependence discussed 
below. 

In the local London theory, at low temperature (typi- 
cally T < T c /3), the temperature dependence of the gap 
can be neglected and A(T) is given by : 



(X(T)-X )/X = ] j^e-^y^ (3) 

where k B is the Boltzmann constant and A(0) the value 
of the smallest superconducting gap at K. Indeed, due 
to this exponential dependence only the smallest gap is 
expected to show up at low temperature. As shown in 
Fig|2] a very poor agreement to the data is obtained if 
the gap parameter is fixed to the weak coupling BCS 
gap of 1.76/u#T c . A much better fit is obtained taking 
A/k B = 7.2 1.2 T c . Note that this gap value does 
not depend on the renormalization constant. The as- 
sumption X(T m i n ) « Ao, where T min is the lowest exper- 
imental temperature, is valid here, since an extrapolation 
of Eq.(3) gives (A(0.5iT) - A )/A < lO" 3 . Since the su- 
perfluid current mainly originates from the Nb sheets for 
H\\c, those measurements hence directly show that this 
small gap is associated with those sheets in NbS2 as well. 

In order to get information on the larger gap, it is nec- 
essary to perform a detailed analysis of the temperature 
dependence at higher temperature. The inset of figE] 
displays the temperature dependence of the-superfluid 



Figure 2: Temperature dependence of the rescaled frequency 
shift Sf/Afo oc AA a b/Ao in NbS 2 at T<T c /3 (red circles). 
Results are fitted with the exponential low temperature ap- 
proximation (see text) to determine the smaller gap A. An 
excellent fit is obtained for A « 1.17 ksT c . Results can't be 
described by an exponential behaviour with the weak coupling 
BCS value A = 1.76 ksTc. Inset: superconducting transition: 
T c =6.05 K is determined by the maximum of variation of the 
frequency shift. 

density p a b(T) up to T c deduced from Eq.(2) for differ- 
ent R/XoAfo values compared to Hf/Hf(0). Note that 
if a small change in R/XoAfo will modify strongly the 
normalized superfluid density, the curvature of Hf(T) is 
utterly independent of Ao- Hf(T) can hence be used to 
determine R/XoAfo with good accuracy and, as shown, a 
good agreement between Hf and the superfluid density is 
obtained taking R/XoAfo = 0.02i^z _1 . Clear deviations 
from the BCS weak coupling theoretical curve (dashed 
line) is visible on the entire temperature range. In par- 
ticular, the BCS values exceed the experimental data at 
low temperature due to the small value of the gap. 

In the multigap model, p is phenomenologically as- 
sumed to be the sum of two independent superfluid den- 
sities 

p = UJpl + (1 - Uj)p2 (4) 
where uo is the weight of band 1 and pi are given by 

The best fit to the data using Eq.(4) and Eq.(5) (solid 
line in Figj3]) is obtained with Ai = 1.05±0.05x ksTc and 
A2 = 1.8 ± 0.1 x ksTc with uj ~ 0.5. Those values are in 
very good agreement with those deduced from tunneling 
spectroscopy measurements which revealed a distribution 
of gaps having two main values at 0.53 meV=6.15 
1 k B T c and 0.98 meV=10.8 1.8 k B T c 15 . Moreover, the 
temperature dependence of the specific heat measured on 
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Figure 3: Normalized super fluid density of NbS2 de- 
duced from Tunnel Diode Oscillator measurements (TDO) 
for RXo/ Afo=0.02Hz~ 1 (closed squares), compared to 
Hf(T)/Hf(0) (open circles) both on sample C8. As shown, 
the weak coupling BCS normalized superfluid density (dashed 
line) is not compatible with the data. A very similar tem- 
perature dependence was previously obtained in NbSe2 (solid 
triangles) and both compounds could be well described by a two 
gap model (thick line, see text for details). Inset: influence of 
the R/XoAfo parameter on the shape of the superfluid den- 
sity, R/XoAfo being equal to respectively 0.01, 0.02 and 0.03 
Hz~ x for diamonds, squares and circles. 



sample C8 can be well described taking Ai ~ 1.05 x ksTc, 
A 2 ~ 2.1 x ksTc and uo ~ 0.6 16 . Even if this two gap 
model well reproduces the temperature dependence of pi 
it is worth noting that this temperature dependence is 
also compatible with a single anisotropic in-plane super- 
conducting gap. The simplest model which respects the 
symmetry of the crystal has an in-plane 6-fold shape, 
A((/>) = Ai + A<2 2 Al (1 + cos(60)) where the extremum 
of the superconducting gap define two characteristics en- 
ergy scales. With this model, we find Ai ~ 0.9(1) x ksT c 
and A 2 ~ 1.9(3) x k B T c . 



IV. COMPARISON BETWEEN NBS 2 AND 
NBSE 2 

We have also reported in Fig |3] the temperature depen- 
dence of deduced from penetration depth measure- 
ments in 2H-NbSe2 12 . Since the critical temperatures 
are slightly different in the two compounds, the data have 
been plotted as a function of T/T c with T c = 6.05 K and 
7.1 K for NbS 2 and NbSe 2 , respectively. In both case, the 
presence of two characteristic gap values are necessary to 
fit the temperature dependence of the in-plane superfluid 
density. We hence show that in both compounds the ex- 
otic superconducting gap structure is related to the Nb 
tubular sheets and that, even if the charge density wave is 



perturbing those sheets in NbSe 2 (but not in NbS 2 which 
does not exhibit any CDW), this CDW does not affect 
the superconducting gap structure. Our results hence 
contradict the ARPES measurements on NbSesri which 
suggested a strong correlation between the gap structure 
and the CDW formation. The presence of an identical 
exotic gap structure in two different compounds (even 
of the same family) is particularly rare which highlights 
the robustness of the underlying mechanism at its ori- 
gin. It's worth noting that there exists a R/XoAfo value 
(~ 0.02i^z _1 , see above) for which the temperature de- 
pendence obtained on both compounds are very similar 
on the whole temperature range. 

Note that this strong similitude in the gap structure is 
not in contradiction with the fact that the vortex cores 
have an isotropic shape in NbS 2 15 , and 6-fold shape in 
NbSe 2 32l3S . Indeed, in this later case, the particular 
shape could be due to an anisotropic Fermi velocity 34 
caused by the presence of the Charge Density Wave. 

Finally, It is interesting to note that even with sim- 
ilar superconducting gaps, it is possible to understand 
the large difference of the anisotropy of the upper crit- 
ical field in the two systems. In the case of NbSe 2 , 
the anisotropy of the upper critical field can be under- 
stood by the anisotropy of the two set of Nb-bands. In- 
deed, the coherence length related to the pancake band 
is larger due to a large renormalised Fermi velocity and 
a small gap 8 . According to band structure calculations, 
the anisotropy of the Fermi velocities (v^* /v c F ) of the two 
tubular sheets are on the order of 2 (Band 17—) and 10 
(Band 18) in NbSe 2 , none of which corresponding to the 
experimental value of the anisotropy of the upper crit- 
ical field (~ 3.2) 23 . However, /i #J 2 (°) = $o/2< 2 6 is 
defined by the band having the smallest ^ value. All 
bands have quite similar in plane Fermi velocities but if 
one assumes that A 2 is associated with band 17 and Ai 
to band 18, on gets f aM7 ~ (Ai/A 2 )f a 6,i8 ~ faMs/1-7 
and H^2 will be defined by band 17. Moreover, for 
/io#J 2 6 (0) = 3>o/27r£ a b£ c , and the very small value of 
c-axis Fermi velocity of band 18, will imply that i^J 2 5 
will be defined by this latter band. One hence gets 

Tj\\ab / -rj\\c ( ab,17 ab, 17 / ab,18 c,18\ w /A / A \2 

H c2 l H c2 = ( V F - V F l V F - V F ) X ( A l/ A 2) 2 ~ 

10/(1. 7) 2 ^ 3 in good agreement with the experimen- 
tal value. Note that this situation strongly differs from 
the one observed in MgB 2 in which the anisotropy of H C 2 
is defined only by the a-band (for T «T C ). As NbSe 2 

and NbS 2 have similar i^J 2 6 values, one can conclude 
that band 18 is not affected by the chemical substitution 
and the increase of the H C 2 anisotropy to ~ 8.1 in NbS 2 
is a hence due to an increase of v^' 17 but this change in 
the Fermi velocity does not affect the superconducting 
gap structure. 

The presence of two characteristics gap values appears 
to be a robust characteristic feature in 2H dichalco- 
genides, but its origin is still open to questions. With 
the present results, the influence of the CDW also has 
to be excluded, underlying the necessity to propose new 
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clues for the origin of the reduced superconducting gap 
observed in both compounds. 

V. CONCLUSION 

To conclude, the temperature dependence of the in- 
plane magnetic penetration depth of single crystals of 
2H-NbS2 has been measured, and compared to previous 
measurements on 2H-NbSe2- The temperature depen- 
dence shows a similar curvature in both compounds evi- 
dencing a similar reduced superconducting gap structure 
in the sheets of the Fermi Surface where the e-ph cou- 



pling constant are the largest. The CDW, present only 
in NbSe2, is not at the origin of the reduced supercon- 
ducting gap. In the light of the discovery of a particular 
acoustic phonon mode observed in lead or Niobium 36 , 
a more accurate knowledge of the phonon dispersion in 
NbS2 could shed insight into the origins of the reduced 
superconducting gap in the dichalcogenides. 
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